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1. Introduction 

The strongly interacting non-Abelian gauge theory is a central key-concept in the high energy 
physics. One of the most important properties is the asymptotic freedom, which indicates a phase 
transition (or crossover) from hadron phase to a deconfined matter -Quark-Gluon Plasma (QGP)- 
in Quantum Chromodynamics (QCD) when temperature {T) and/or quark chemical potential (/x) 
exceed a typical energy scale of QCD (Aqcd)- To investigate the QGP and this phase transition is 
a physics goal in the LHC-ALICE experiment, and indispensable to understand the astrophysical 
systems, such as the early universe and the neutron stars. 

The phase transition (or crossover) from the hadron to the QGP can be described in terms 
of two kinds of (approximate) symmetries, the SU(A'y ) x SU(A'y ) chiral symmetry and the center 
Zat^ symmetries {Nc and Nf represents the number of colors and flavors). The finite value of the 
Polyakov loop indicates a sort of deconfinement in the sense that an excitation energy of heavy 
quark is finite. The finite value of the chiral condensate at vanishing current quark mass indicates 
the hadron phase in the sense that hadrons are massive (the pion is massless) and relevant degrees 
of freedom. The Polyakov loop and the chiral condensate breaks the Zat, and the chiral symmetries, 
respectively, and are typical observables in the lattice QCD. The lattice QCD based investigation 
for the phase structure of their symmetries i.e. the QCD Phase Diagram in T — /i plane is one of 
the most important subjects in the high energy physics. 

The Lattice QCD Monte-Carlo simulations (LQCD-MC) indicate that chiral and deconfine- 
ment (ZjvJ crossovers simultaneously happen around, =~ 160— 190 MeV at = [[T]]. The co- 
incidence of "Zat, and chiral" is the non-trivial observation, and the relation between their crossovers 
is a modem interest on the QCD phase diagram. The important sampling method in LQCD-MCs 
suffers from the sign problem of the quark determinant with finite }X, and the LQCD-MC does not 
work well at finite jx. It is required to shed light on the large pL region of the QCD phase diagram 
beyond the sign problem. 

In this proceedings, we investigate QCD phase diagram in the whole range oi T — jx plane 
by using the Strong Coupling Expansion in Lattice QCD (SC-LQCD) with one species of (un- 
rooted) staggered fermion. The analytic investigation based on SC-LQCD can be informative for 
finite jX, and there are many recent developments ||3|-^]. We take account of next-to-leading order 
(NLO, ^(1 /g^)) terms in the strong coupling expansion. In addition, we consider the leading pure 
gluonic contributions for Polyakov loops. We concentrate on the leading order of the large dimen- 



sional {l/d) expansion [ |10| ] for simplicity. In these setups, we focus on the relation between the 



U^{\) chiral and Zj^^, deconfinement dynamics in the phase diagram. 



2. Effective Potential 

We start from the lattice QCD partition function with one species of staggered fermion {%) 
with a quark mass (mo). Gluons are represented by hnk variables (?7v,x)> 



^r]v.xXxUv,xXx+v-nvl{h.c.) Up + Up 

L— L"'o{xx)x+2^ — -r^ 

VJC ^ X p 6 



(2.1) 
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where, [/p=^v,.r = tv,[U^,^Uy^+^ul^^^ul,\, and (t]o,.v, T];,.,) = [e^ ,[-iy^+-+^'-^). By using a 75- 
related factor Ex = (— l)^o+ '+^'', a staggered chiral transformation is defined as Xx — ^ e'^^^Xx [@]- 
The lattice action is invariant under this U^{\) chiral transformation in the chiral limit /mq — 0. In 
the following, we use the lattice unit a = \, and we investigate the color S\]{Nc = 3) case in 3 + 1 
dimensions (J = 3). 

We perform strong coupling expansions, and take account of the next-to-leading order (NLO) 
effects. To manipulate temperature effects, we firstly evaluate the spatial gluon degrees of free- 
dom by utilizing the group integral formula, JdUjx Ujx'^Pj^^ = 5^8^. /Nc, which leads to sums of 
color singlet (hadronic) composites. The spatial link integral leaves an isotropic sum over spatial 
directions, and the energy per bond is assumed to be proportional to l/din order to keep the action 



finite at large spatial dimension d |[10p. In this normalization scheme, the quark field scales as 
J- 1/4 We concentrate on the leading order \/d terms, which correspond to the minimum quark 
number configurations for a given plaquette geometry. Within these approximations, we obtain the 
following hadronic composites and effective couplings, 

{Mx, Vx+{^), Vx-m = iXxXx, e^XxUo,xXx+6^ e-^Xx+oKx^-r) , (2.3) 

where the lattice coupling j3 = 2Nc/g^ is the unique parameter in the current theory. The left 
panel of Fig. |l| provides a schematic description of composites. In addition to them, we extract the 
leading pure gluonic effect to Polyakov loops (L^ x = nT^o,A > the right hand side of Fig. |I|), which 
is necessary to investigate a deconfinement dynamics. These setups give us the simplest framework 
to investigate both Uj^{\) chiral and Za?^ deconfinement transitions (or crossovers) in the SC-LQCD 
with finite coupling effects. Recent our study indicates that the next-to-next-to leading order effects 
with quarks has just a small correction to the phase diagram, while there is discrepancy between SC- 
LQCD and the Lattice QCD Monte-Carlo simulations (LQCD-MC) for the jS = 2Nc/g^ evolutions 
of the critical temperature at zero chemical potential, (Tc,^=o{^)) [||]. Due to the pure gluonic 
Polyakov loops, the Tc,^=o{p) is found to become closer to the lattice Monte-Carlo (LQCD-MC) 
simulations around j3=4.0at/i=0[^. We anticipate that the pure gluonic Polyakov loops could 
have an essential contribution at "finite /x", which will be investigated here. 

To integrate the quark degrees of freedom iX,X)^ we introduce four kinds of mean fields |^], 
which are summarized in Table The mean field a represents the chiral condensate, which leads 
to the dynamical mass shift, niQ — = niQ + {d / {2Nc))o in the strong coupling limit It is remark- 
able that the NLO effects are expressed as modifications of quark mass and chemical potential, 
(m^,/x) — )• {mq,jx) with the quark wave function renormalization factor, yJZ+Z^ , and those are 
governed by the unique lattice parameter j8 [^, 

^ '"o + + 2As(p.v) a ) , At = M - log W , (2.4) 



where, Z± = 1 + j8 ((Pt ± Wt) ■ 

This NLO formulation which is invented in our previous work allows us to evaluate the 
remnant degrees of freedom, quarks ix,X) ^^d the temporal link variable Uo, in the same way to 
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the strong coupling limit case [11, 12]. It is convenient to take a static and diagonalized gauge 
(called the Polyakov gauge) for temporal link variables with respect for the periodicity [|ll|], ?7o,x = 



diag{ 



} . In this gauge, the Polyakov loop reduces to a simple expression, Lp x 



: Nclp,x, and the Haar measure is also expressed in terms of the Polyakov loops. 



Haar(^p,X) ^p,x) 



1-61 



(2.5) 
(2.6) 



The Haar measure effects are responsible for the deconfinement dynamics. The quark determinant 
includes the quark propagation wrapping the temporal direction with Uq, and leads to the quark 
driven Polyakov loops, which is characterized by, 



1 +A^^(/^,^^-(£,-A)/r + /^^^^-2(£,-A)/r) _^^-NAE,-fi)/T 
1 +iV^.(/^_^^-{£,+A)/r + ;^^^^g-2(£,+/i)/r) +^-NAE,+fi)/T 



(2.7) 
(2.8) 



Here, Eq{mq{o)) = sinh^' [m^(a)] corresponds to the quark excitation energy. The Polyakov loop 
Ip couples to a Boltzmann factor ^^(^9^'')/^, and determines quarks thermal excitations. Note that 
the Boltzmann factors is a function of the chiral condensate a, therefore the chiral and deconfine- 
ment dynamics couples to each other via Eq. (IH) and (2^). This effect has been investigated 



in the strong coupling limit case []12p, and now we have the NLO corrections via Eq. ([2.4|). In- 
stead of integrating out the temporal link variables, we replace the Polyakov loop with its constant 
mean-field value, {lp.xjp,x) (Ipjp)- Finally, we obtain the effective potential as a function of the 
auxiliary fields = (a, (p^,*, COtJpJp), temperature T, and quark chemical potential /i. 



- N,Eq - T (log %{T,n) + log^,-(r, At )) 



2dTN^ 



l/T 



Iplp - riog ^Haar(^p, ^p) • 



(2.9) 



The equilibrium is determined by imposing stationary conditions on the effective potential, 5^eff/^*J* - 
0, which leads to the relations summarized in the third column of Table |l|. 



C + 60 

V+ K 



1/,'/' ii 
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Figure 1: Left: The strong coupling limit and effects to appear within the leading order of l/c/ expan- 
sion. Open circles, filled circles, and arrows show t/y, respectively. Right: Schematic figure of the 
pure gluonic Polyakov loop in the strong coupling expansion. 
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Table 1: The auxiliary fields and their stationary values. Here, = Nc — ^/Z+ZZfhq + jixCO^. 



Aux. Fields 


Mean Fields 


Stationary Values 


a 


i-M) 


-dr,/dWZ+Z-'ng) 


(ps 


(MM) 






-((y+-y-)/2) 


2(po/il + y/\+4l5r(po) 


(Or 


-((y+ + v-)/2) 





3. Results 

We will concentrate on the j8 = 4.0 results in the followings. In the left panel of Fig. |^ we 
show the phase diagram in the chiral limit. In the low T and large pL region, we find the critical 
end point (CEP), and observe a partially restored (PCR) matter around the CEP. The details of PCR 
have been explained in our previous paper [^. In this proceedings, we concentrate on the relation 
between the chiral and deconfinement crossovers. In the right panel of the Fig. ^ the equilibrium 
value of the Polyakov loop Ip is depicted as a function of Ta and jia in the same condition as the left 
panel. We find the rapid variation of Ip around the chiral transition line. Thus the deconfinement 
and chiral dynamics are strongly correlated except zero temperature cases. 

Let us go along with the jia = 0.5 (dashed blue) line in the right panel of the Fig. ^ The 
derivative of the Polyakov loop in terms of T on this line is shown in the left panel of Fig. |3[ In 
the chiral limit (mo — )• 0, left panel), we find two peaks "P" and "Q", and the former corresponds to 
the steepest point on the pLa = 0.5 line in the right panel of Fig. ^ (Here, "P" and "Q" in the Fig. || 
correspond to those in the Fig. |2|). Near the first peak "P", the Polyakov loop is force to change 
its value to minimize the effective potential due to the drastic variation of the chiral condensate. 
In this meaning, the peak "P" can be interpreted as the Chiral Induced Deconfinement Crossover. 
When the current quark mass mo becomes large, the Zn^ center symmetry becomes dominant. The 
right panel in the Fig. ^ corresponds to such a situation (moa = 1.0). We find that the second peak 
"Q" grows up, and the "P" disappears. This indicates that the second peak is the Z/v^ Induced 
Deconfinement Crossover. 

It is interesting that the Za?^ nature still survives in the chiral limit, and its pL dependence is 
almost negligible. Hence the chiral and Z^v^ dynamics start to separate with increasing /i. This is 
the mechanism of the two sequential deconfinement crossovers. We note that the coincidence of 
the chiral and Zj^^ at "/x = 0" is still mysterious, and to be investigated in future. 

4. Summary 

We have investigated the chiral and deconfinement crossovers at finite temperature {T) and 
chemical potential {p.) based on the strong coupling expansion in the lattice QCD with one species 
of the staggered fermion. We have considered the next-to-leading order (NLO) effects in the strong 
coupling expansion with the leading pure gluonic contributions to the Polyakov loops. We have 
concentrate on the leading order of the \/d expansion, and derived the analytic expression of 
the effective potential in the framework of the mean-field approximation. We have replaced the 
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Figure 2: Left: The phase diagram in the chiral limit at j3 = 4.0. Right: The equilibrium value of the 
Polyakov loop as a function of T and /i in the chiral limit at )3 = 4.0. 




T 



Figure 3: The derivative of the Polyakov loop in terms of T at /i = 0.5 in the lattice unit. The left and right 
panels show the chiral limit and heavy mass (moo = 1 .0) cases, respectively. In the left panel, the vertical 
dashed (blue) Une represents chiral transition temperature. 



Polyakov loop with its constant mean fields value without integrating the temporal link variable, 
and the deconfinement dynamics is introduced via the Haar measure of temporal link integral. 

At a fixed large /i, we have observed sequential two deconfinement crossovers i.e. the rapid 
variations of the Polyakov loops as T increases. One of them results from the strong correlation be- 
tween the chiral and deconfinement dynamics. In other word, the Polyakov loop is force to change 
its value to minimize the effective potential due to the drastic variation of the chiral condensate. 
In this meaning, we have found the Chiral Induced Deconfinement Crossover. When the current 
quark mass mo becomes large, this crossover disappears, and the another crossover grows up. Thus 
we have also observed the Zj^^ Induced Deconfinement Crossover. In the chiral limit, the former 
becomes stronger, and the latter becomes weaker. The point is that the latter still exists in the chiral 
limit. Hence two crossovers co-exists there. 

There are several studies to be performed in future. Firstly, it has been known that the tem- 
poral link integral can be analytically performed, and then the fluctuation effects appear as the 
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combination of the modified Bessel function and the difference of Polyakov loops in the effective 
potential It should be confirmed that the obtained results in this work is independent of approx- 
imation schemes. Secondly, the chiral and Polyakov loop susceptibilities would be informative to 
investigate the relation between the chiral and deconfinement dynamics. Thirdly, the exact eval- 
uation in each order of the strong coupling expansion is required beyond the \/d expansion and 
the mean-field approximation. This could be achieved by extending the Monomer-Dimer-Polymer 



formulation |n, 13] to include the NLO effects. And finally, it would be interesting to introduce the 



advanced formulation of the pure gluonic effects for the Polyakov loops which is recently investi- 



gated in Refs. [14|. 
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